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Abstract 

We consider the quantum (trajectories) filtering equation for the case 
when the system is driven by Bose field inputs prepared in an arbitrary 
non-zero mean Gaussian state. The a posteriori evolution of the system is 
conditioned by the results of a single or double homodyne measurements. 
The system interacting with the Bose field is a single cavity mode taken 
initially in a Gaussian state. We show explicit solutions using the method 
of characteristic functions to the filtering equations exploiting the linear 
Gaussian nature of the problem. 


1 Introduction 

The aim of this paper is to write explicite solutions to equations for the quan¬ 
tum trajectories for a linear system (cavity mode) driven by inputs in a general 
Gaussian state. The quantum filtering theory developed by V. P. Belavkin [T]- 
|T] describes the dynamics of a dissipative Markovian quantum system under an 
indirect and continuous in time observation. The model is a quantum analog of 
the classical filtering theory [5]-[10] with the framework of the quantum stochas¬ 
tic Ito calculus (QSC) [TT| being necessary to capture the physical noise. The 
Bose field interacting with the system here plays the role of measuring appara¬ 
tus. The evolution of the system conditioned by the results of the measurement 
is given by the filtering equation, and the conditional density matrix depending 
on the random results of the measurement can be viewed as a stochastic process. 
The theory of quantum filtering, or quantum trajectories as it is known in the 
physics literature, has been developed and is used nowadays in a technique of 
quantum measurement and control |12j - m- 

The Belavkin filtering equation, also known as the stochastic master equa¬ 
tion, was derived for the Bose field in a pure Gaussian state such as a vacuum 
|3| , squeezed vacuum (TS] , and coherent state [281 . The rigorous derivations 

of the quantum filter in more general case when the systems is coupled to the 
multiple fields in arbitrary zero-mean jointly Gaussian state was given in m- 
The filtering equation for the case when the output field mixed with the field in 
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a Gaussian state one can find in m- In the latter we presented the result for ar¬ 
bitrary zero-mean and nonzero-mean Gaussian state of the Bose field including 
vacuum, squeezed vacuum, coherent, thermal, squeezed thermal and squeezed 
pure state. Recently the quantum filtering theory has been also expanded to 
the non-classical states of the Bose field such as a single photon states [31j and a 
superposition of continuous-mode coherent states [33] . We mention that recent 
experiments have shown that single photon homodyning may become a realistic 
procedure |34| . 

In this paper we describe the a posteriori evolution of a harmonic oscillator 
(an optical cavity mode) driven by the Bose field in a Gaussian state for two 
measurement schemes: a single homodyne measurement and a double homodyne 
measurement. Let us remind that in the first of these schemes we measure only 
one quadrature of the system while the second one allows us inaccurate joint 
measurement of two quadratures of the system. In the paper we provide the 
reader with the differential equations for the a posteriori mean values of the 
system’s operators, show the condition for preservation of the purity of the 
system’s state and assuming that a cavity mode is initially in a Gaussian state 
we give an analytical solutions to the quantum filter and compare them with 
the results for a priori dynamics. 


2 Quantum stochastic processes 

Let us consider two independent annihilation processes Bj{-) {j = 1,2) and 
related to them two independent creation processes Bj (•)* (j = 1,2) satisfying 
the commutation relations m-m 

[B, (t) , Bk (s)] = [B,ity, Bk{s)*] = 0 , [B, it) , Bk (s)*] = S,k (t A s) , 

where t A s = min(t, s). The processes Bj (•) (j = 1,2), Bj (•)* (j = 1,2) are 
quantum analogues of classical Wiener process. 

We analyze the case when the baths are in Gaussian states with the quantum 
mean values 

E [Bj (t) Bk (s)] = 6jkmj {t A s) , 

E [Bj (t)* Bk (s)] = Sjkrij (tAs) , 

E [Bj (t) Bk (s)*] = 6jk{nj + 1) [t A s) , 

E[Bjit)]=Pjt (j = l,2). 

The ltd multiplication tables are given then by 


X 

dBjit) 

dBjit)* 

dBjit) 

rrijdt 

irij + l)dt 

dBjit)* 

rijdt 

m*dt 


where Uj > 0, \mj\^ < rij {rij + 1). All the other products, including dBi{t), 
dB 2 {t)*, etc., vanish. 
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When j3j = 0, nij = 0, rij = ^^ where ojj is the carrier 

frequency of the j field, we get a thermal bath. In general, the processes Bk 
may be represented by an Araki-Woods double-Fock space construction, as done 
explicitly in [35]. However, in the degenerate case where we have = 

Uj {rij + 1) one can use for each k the representation = •^/n + lAk + \/nA^ 
where Ak is the annihilation process on a single copy of Fock space, with the 
choice of Fock vacuum. For a finite number of states, the latter case corresponds 
to Gaussian pure states (vacuum, squeezed vacuum or squeezed coherent state) 

[Ml Eg. 


3 Interaction between observed system and mea¬ 
suring apparatus 

We assume that only the first input, Bi{t), interacts with the quantum system 
in interest (system S) and the evolution of the compound system is described by 
the unitary operator, C/(t), which satisfies the quantum stochastic differential 
equation (QSDE) of the form [11] [STj 

dU{t) = {LdBi{t)*-L*dBi{t)-{K + iH)dt}U{t), U{0) = I, 

where 

^ ^ (ni + 1)L*L + iniLL* - , 

H is the Hamiltonian of S, and L is a coupling operator. 

The above equation is written in an interaction picture rotating at the fre¬ 
quency wi and H is the only part of the Hamiltonian not removed by this 
procedure. 

For any observable X of the system S in the Heisenberg picture 

jt{x) = u{tY{x®i)u{t), 

we obtain the QSDE [3S] 

djYX) = {-ijY[X.H]) 

+ 2 A [X, L] + [L*,X]L) 

+ ^nut {L[X,L*] + [L,X] L*) 

-^mljt{L[X,L] + [L,X]L) 

-imi jt {L* [A, L*] + [L\X] L*)} dt 

+3t{[X,L]) dBYty +jt{[L\X]) dBYt), (2) 

which one can check it by making use of the rules of QSC and the multiplication 
table ([T|). 
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Let us remind that according to the interpretation of the model given by 
Gardiner and Collet [40j the operator Bi{t)* refer to the Bose field be¬ 

fore the interaction with the system S, while the field after the interaction is 
described by the output processes 

= U{ty {I ® Bi{t))U{t), 

B°'^\t)* = Uit)* {I®Bi{t)*)U{t). 

Thus the output field is the field in the Heisenberg picture and we have [22] 
Br\t) = [ {jsiL)ds + dBys)) , 

B°yyty = [ Us{L*)ds + dBys)*). 

Jo 

The output field conveys information about the system S. Note that the incre¬ 
ments dBy^^(t), dBy^^yy satisfy the multiplication table of the form (|T]). Thus 
the output field remain Bose free field. 

4 Posterior dynamics for a single homodyne ob¬ 
servation 

In this section we shall describe an evolution of a single-mode field in a cavity 
conditioned by the results of a continuous in time measurement of an optical 
quadrature when the cavity is driven by a vacuum, coherent, thermal or squeezed 
input. 

We take the hamiltonian H = hda*a, where 5 = Wc — wi, Wc is the frequency 
of the cavity mode, and a stands for annihilation operator, and the coupling 
operator L = y/Jla with ^ G K+. 

Let us consider the continuous measurement of the output process 

Y{t) = (3) 

where 9 G [0,27r). Observation of Y{t) can be interpreted as the indirect and 
imperfect measurement of the cavity mode observable -I- 

A measurement of phase-dependent quantities is possible by making use of 
the homodyne/heterodyne detection scheme where a signal escaping from a cav¬ 
ity is superposed with an auxiliary field (a local oscillator) using a beam splitter 
I21I3ZI- In practise, the measurement of Y (t) is realizable by an implementa¬ 
tion of the homodyne scheme. We get the homodyne measurement if the carrier 
frequency of the input field is equel to the carrier frequency of the local oscil¬ 
lator. When this condition is not satisfied, we have the case of the heterodyne 
measurement. 
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The density matrix of the cavity mode under the continuous observation of 
the process V (t) satisfies the stochastic filtering equantion of the form [31] 

dfft = ^-i[5a*a,Qt] +y/Jl[a,Qt] j3l+y/Jl[Qt,a*] Pi 

+ ^ + 1 ) {[agt,a*] + [a,gta*]) 

+ {[a*gt,a] + [a*,gta]) 

-^rnl {[agt,a] + [a, gto]) 

-^mi ([a*£>t,a*] + [a*, gta*]) | dt 

+ ^{e*®oet + e~''^gta* - (e*®(a)t + e“*®(a*)t) gt 
+ (e‘®ni + [a, gt] 

+ {e~"^ni + e^^mi) [gt,a*]^dY {t) , (4) 

where 

dY{t) = dY{t)-^{e^^{a)t + e-^^{a*)t)dt 
-(e*%+e-*®/3r)dt, 

K = 1 + 2ni + and {.)t = Tr(pt.) is the a posteriori mean 

value of an operator. The above equation describes the condition evolution 
of the cavity mode with respect to the commutative von Neumann algebra 
y* = {Y (s) |s < t}. Note that the mean value of dY(t) is zero, {dY(t))'^ = ndt 
and the process Y{t)/y/K = J*dY{s)/^/K is izometric to the standard Wiener 
process. 

One can check that if \mip = ni(ni + 1), Eq. (|4|) transforms pure states 
into pure states and it is equivalent to the stochastic filtering equation 

d\'4>t) = { — id a* a — ^ (ni + 1) a* a — ^niaa* 

+ ^mi(a*)^ + + y/JiPla - y/JiPia* 

+ — - {a*)t [a (e*®ni + + e*®) 

K 

-a* (e*®TOi + e“*®ni)] - ■^{a)t{a*)t}\P’t) dt 

+ To (e'^ni + + e*®) 

K 

- a* (e*®mi + e“*®ni) - e*®(a)t] \pJt) dY{t ), 

for the a posteriori wave function \^pt)■ Here (.)t = {tptliPtpt) and initial condi¬ 
tion states as \'ipt=o) = IV'o)- Thus we obtain preservation of state purity when 
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the field driving the system is pure. The validity of the above equation one can 
check by differentiating gt = |'0t)('*/'*I- 

To solve Eq. we will use the normal ordered characteristic function 

associated with the state: 

0(r,e) = Tr(e-*«*>e-*«“*) (5) 

with the inverse given by 

Note that p is related to the P-function in quantum optics 127]: here the state 
is represented in terms of coherent states as p = f |a)(Q;| P{a)(Pa and P(a) = 

Making use of the transformation ([S]) , one can convert Eq. (|2|) into the form 

+y/JdldiiC + y/JdPliS, - pm 

(e“*®TOi+e*®ni)-e*®(a)t-e“*®(a*)t] g{C,i;t)dY{t),{6) 

where = d/d^, = d/d^*. 

We assume that the cavity mode is initially in the Gaussian state 


HC,^,t = o) 


exp 


-i(rao+eaS) 




(7) 


where OQ) Co G C, uq > 0. One can easily check that 

Tr[apt=o] = Qfo , Tr[a^pt=o] = Co + : Tr[a*apt=o] = t'o + |ao|^ • 


We shall prove that the solution to Eq. ([21) corresponding to the initial state 
© has the Gaussian form 


p(r,C;i)=exp 


-*(r(a)t+C(a*)t) 




( 8 ) 


where {a)t = Tr[opt]. 

To this end, it is convenient to rewrite Eq. in terms of the stochastic 
function 

= -inp(r,C;t). (9) 

Using the Ito formula 


dl = ——dp + 


m- 


{dp) 
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and the tables o, we obtain 


-y/JiPiiC - y/Jj-Pm + Mni ICI^ 

2 2 J 

+ A I + e-^^id^ I - iC (e*Vi + e-^^m) 

ZAC 

-i^ + e*®ni) + e'^{a)t + e“*®(a*)t] ^ dt 

+ I + e~^^id^ I — iC (e*®TOi + e“*^ni) 

AC 

+ e*®ni) + e''^{a)t + e“*®(a*)t] (iy(t). (10) 

And finally inserting (jO]) into Eq. m and equating expressions of the same 
power of parameters ^ and provides us with the consistent set of the differ¬ 
ential equations: 

d{a)t = - (iS + ^ {a)tdt - ^/JiPidt 

_ 1 _V^ [e*® {({t) - mi) -I- - ni)] dY{t), (11) 

AC 

C(t) = -2i6C{t) - n iC{t) - mi) 

-- [e*® (C(0 - mi) + {v{t) - ni)]^ , (12) 

AC 

and 

z>(t) = -/X (i/(t) - ni) - — |e*® (C(t) - mi)-I-(z/(t) - ni)|^ (13) 

with the initial conditions: {a)t=o = oq, C(^ = 0) = ()o, and v(t = 0) = vq] 
which completes the proof. 

One can see from Eqs. (HIHISl) that the coherent state of the cavity mode 
is preserved under the observation of Y (t) only when the input field is taken 
in a vacuum or coherent state. If mi ^ 0, then the a posteriori state becomes 
squeezed. In the resonance case when d = 0, we obtain the stationary asymptotic 
solutions of the form: (a)oo = ; C(oo) = mi, and v{po) = ni. So the system 

approaches the Gaussian state with the parameter of squeezing defined by the 
state of the input field. When mi = 0 then for any value of 5 we have C(oo) = 0 
and v{oo) = ni. 

We may draw from Eqs. (HIHSl) also the conclusion that the a posteriori 
mean values of optical quadratures X = a-I-a* and P = (a — a*) /i of the cavity, 
in general, depend on the measured noise, whereas their dispersions, given by 

(AX(t))2 = 1 + 2iy{t) + 2ReC{t), (14) 
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(AP(t))2 = 1 + 2iy{t) - 2Re({t), 


(15) 


remain non-random. 

Let us note that Eqs. m and (1131) can be written as matrix Riccati differ¬ 
ential equation of the form 


where 


Z{t) = Z{t)TZ{t) + RZ{t) + Z{t)R*+W, 
Z{t) = 


i/(t) C(t) 

at)* I'it) 


T = - 




1 

„2ie 


-2ie 

1 


R = 


W = 


ni +n\— |toi|^ 


{\mi\ — nl) 

with the initial condition Z(t = 0) = 


mi + e ^“(|mi| — nf) 
ni -I- nf — |toi 


(16) 


■ -iS 

0 

_L ^ 

-^+ni+ 

TOi -I- e 

0 

iS 

H— 

TO* -1- 

-f -1- ni -1- e“^*®TOi 


t'o Co 
Co* ^0 

Eq. m can be solved by using the matrix fraction decomposition technique 
|4T] which consists in representing matrix Z{t) in the form 

z{t) = zi{t)z^at). 


The initial values of Zi{t) and Z 2 {t) are constrained by the initial value of Z{t). 
One can check that the differential equations for the numerator and denominator 
matrices are linear and they can be read as follows 


■ Zi{t) ■ 


. z^it) . 



W 

-R* 


Ziit) 

Z2{t) 


The solution to the foregoing can be written as the product 


■ Zi{t) ■ 

- pTt 

■ Z(0) ■ 

. Z2{t) 

— c- 

1 


(17) 


where T = 
matrix of size 


R 

-T 

2 . 


W 

-R* 


is the time-invariant 4-by-4 matrix and 1 is identity 


5 Posterior dynamics for a double homodyne ob¬ 
servation 

Joint perfect measurement of two optical quadratures of a quantum system is 
not possible, but one can analyze their simultaneous imperfect measurement. 


























Using a double homodyne scheme one can mix the field escaping from the cavity 
(the output field) with an auxiliary field and measure simultaneously the two 
signals: 

Ui(t) = — ((5i(t) + (32(0)) 
l2(t) = ^(Pl(t)-P2(t)) , 

where Qiit) = Pi{t) = /i are quadra¬ 

tures of the field escaping from the cavity and <32(i) = ^2(0 + B 2 {t)*, P 2 {t) = 
(i?2(i) — ji describe the squeezed noise field do not interacting with S. 

The observed processes Yiit) and Y 2 {t) satisfy the condition [Yi(t), 12(3)] = 0 
for all t, s > 0 and their increments 

dYi (t) = {dBi (t) + y/Jljt (a) dt + dB 2 {t) -|- H.c.} , 

v2 

dY 2 (t) = —^ {dBi (t) + y/JI jt (a) dt - dB 2 {t) - R.c.} , 

V2i 

have non-trivial correlation expressed through the table 


X 

dYi 

dY2 

dYi 

(1 -1- ni -1- n2 + m'l -|- m^) dt 

(m'l' — m 2 ) dt 

dY2 

(m'l' — m2) dt 

(1 -1- ni -1- n2 — m'l — m^) dt 


where mi = m[ + im”, m2 = m^ -|- im'^ are the decompositions of the squeezing 
parameters into their real and imaginary parts. 

The filtering equation for the a posteriori state, gt, takes in this case the 
form [3T] 

dgt = i[Sa*a, gt] + y/JI[a,Qt] (31 + ^/JI[gt,a*] Pi 
+ ^ (^1 + 1) i[agt,a*] + [a, gta*]) 

+ {[a*Qt,a] Y [a*, gta]) 

-^mi {[agt,a] + [a, gto]) 

-^mi ([a*£»t,a*] -I- [a*, gta*]) | dt 
+ (1 + ni -f U2 - rn[ - m'^) Ji{gt) 

+i (m" - m2) J2{gt)'^dYi (t) 

2 ^{ (1 + ’^1 + ^^2 + m'l -I- m2) J 2 {gt) 

-i (m" - m2) Ji(pt)|'^^2 (t) , ( 18 ) 
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where 


Ji{Qt) = agt + a* Qt - {a)tgt - {a*)tgt 

+ {ni + ml) [a, gt] + {ni + toi + 1) [gt, a*], 


and 


J2{gt) = agt - a*gt - (a)tgt + {a*)tgt 

+ {ni - ml) [a, gt] - {ni - mi + 1) [p*, a*], 

dYi{t) = dYi{t) - {{a)t + {a*)t)dt - V 2 {RePi+Re/32) dt, 

dY 2 {t) = dY 2 {t) + ij^i{{a)t - {a*)t)dt - V2{lmPi-lmP 2 ) dt, 
A = (1 + m + 712 )^ — {m'l + m' 2 )^ — {m'l — m' 2 )^ . 


(19) 


Here the condition {pt+dt) = Pt+dt if {pt) = Pt is satisfied if and only if 
= 7ii(ni + 1) and |m 2 p = 712(712 + 1). It means that the conditional state 
of the cavity mode preserves its purity if and only if both electromagnetic fields 
used by us for observing the system are in pure states. For instance, if the first 
field is in a coherent state (/3i ^ 0, mi = 0, tii = 0) and the cavity is initially 
in a pure state, then Eq. dT51) is equivalent to the stochastic filtering equation 

dIV't) = [iS + ^ a*a +pa{a*)t-^{a)t{a*)t 


+\/p{Pla- Pia*) \\'4)t)dt 


y2A 


(1 + 772 — m 2 ) dYi{t) 


-i (1 + 712 + m 2 ) dY 2 {t) \ {a - {a)t) It/’*) 


for the a posteriori wavefunction It/’*)- 

We shall prove that Eq. HID transforms Gaussian states into Gaussian 
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states. Using ([5]), we obtain from (fT51) the stochastic equation 


+y/JIl3iiC + - g-ni 1^1^ 



+i^iml + ni) - {a)t - {a*)t] giC,^',*) 
X [(1 + ni + n2 — m'l — m'2)dYi{t) 

— (to" — m!2)dY2{t)] 


H—- idi + iCi'mi - ^i) 
V2Ai 

- m) - (a)t + {a*)t] g{C,C,t) 
X [(1 + ni + 712 + m'l + m'2)dY2{t) 

— {m'l — m'l)dYi{t)\. 


( 20 ) 


In this case the equation for the function l{t) has the form 


-y/JdldiiC - + g-ni |C|^ 



-iC{ml + ni) + {a)t + {a*)t] 

X [(1 + rzi + 712 — m'l — m'2)dYi{t) 

— {m'l — m'l)dY 2 {t)] 

+ ~ 'iC{mi — Til) 

+i^{m*i - m) + {a)t - (a*)t] e(r,CU) 

X [(1 + Til + 712 + m'l + m'2)dY2{t) 

-{m'l - m'l)dYi{t)] . (21) 


Now inserting 


i{c,^-,t) = -t{e{a)t+^{a*)t)-l{c^c{t)+ec{tr)-\eHt) 


into Eq. (^11) and evaluating {dl)^ we end up with the consistent set of differ- 
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ential equations: 


d{a)t 



{a)tdt - y/jiPidt 

[(1 + ni + n2 — m[ — m2)dYi{t) — (to" — 7712)^^2 (t)] 

X (C — TOi + 7/ — m) 

—i[(l + ui + n2 + m'l + m'2)dY2{t) — (to" — m2)dYi{t)] 
x((C- wi)-(7/-ni))|, (22) 


C(t) = - 2 i 5 C.{t) + 

+ -£| - 2(1 + ni + n2)(C - TOi)( 7^ - ni) 

+ (toi + TO 2 )[(C - mif + {v- mf] 

- TO2)[(C - mif -{ly- ni)^]| , (23) 


7>(t) = -n{iy{t)-ni) 

+ £{ - (1 + ni + n2)[(C - TOi)(C* - ml) + {v - rq)^] 

+ {m[ + m 2 ){v — ui)(C + C* “ 2ml) 

- m2){v - ni){( - (* - 2iTO")| , (24) 


with the initial conditions {a)t=o = ctQ, = 0) = ^ 0 i and v{t = 0) = vq. 

It is easy to see that if d = 0 we have the asymptotic solutions: (a )00 = , 

C(oo) = TOi, and 1 /( 00 ) = ni. 

For the coherent input Eqs. (1^51) and form the homogeneous matrix 
Riccati differential equation 


where 


Z{t) = Z{t)TZ{t) + RZ{t) + Z{t)R* 


Z{t) = 


v{t) C{t) 

at)* at) 


T = 




A 

R = 


- (1 + 772 ) 


irin 


m2 -(1 + ^ 2 ) 

0 


-^-zS 

0 


-^+z5 


and initially Z{t = 0) 


Vo Co 
Co* VO 


(25) 
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Using (HZD we obtain the exact expressions of the parameters: 


at) 


g-(/x+ 2 ii 5 )t 

m 


Co + 


fim*2 - ICol^) 
A (/i — 2id) 




( 26 ) 


where 

m 


'{t) = 


-fit 


m 


1^0 


(1 + 712) - ICol^) 


A 


( 1 - 




(27) 


(l + n2)^k^-|Co|1 
1 +-- - 


(1 _ (1 - 


A2 ^ ^ ' A 

(, „_(y_2z5)A MW 2 C 0 


A (p — 2i5) 

|2 t,,2 




( 


1 — e 


-{fi-\-2i5)t 


A (p + 2ii5) 

(1 _ »-<»-»«.) (1 . 


Let us notice that for the coherent and vacuum inputs (toi = m 2 = 0 and 
rii = 712 = 0) and the cavity mode being initially in a coherent state, we get 
C(t) = 0 and i'{t) = 0 for any t. Thus, in that case the amplitude {a)t is 
independent of the measured noise. The same conclusion we obtain for the 
conditional evolution under a single homodyne observation. 

For the large t the parameters ({t) and v{t) go to zero and the systems 
approaches asymptotically coherent state. Therefore after some transient time 
any information about the initial conditions is lost. 


6 A priori evolution 

In this paragraph we would like to compare the a posteriori dynamics of the 
cavity mode with the corresponding a priori dynamics. A non-selective dynam¬ 
ics of the system we obtain by averaging both sides of Eq. or respectively 
m over the outputs. It gives us the equation 

&t = -i[6a*a,at] +aj^[a,at] PI + aJtat,a*] I3i 

+ ^ (ni -I- 1) ([acTt, a*] + [a, ata*]) 

+ ([a*crt, a] -I- [a*, ata]) 

-^mj {[aat,a] + [a, ata]) 

-^mi {[a*at, a*] + [a*,crta*]) , (29) 


for the a priori state cr*. 
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Let us remind that this equation does not preserve in general the purity of 
the state. 

By transformation ([5]), we get from (1^91) the following equation 


+,/JiPiiC + - pn-i ICI^ 


2 ’ 2 , 

The solution to Eq. (1501) for the initial Gaussian state ([7]) has the form 

1 


CT(r,C;i) = exp 


-*(r(«)*+e(a)t) - - {c\{t)+eatr) - 


(30) 


(31) 


The parameters of the state cr* appearing in (1311) satisfy the set of ordinary 
differential equations 


{a)t = - + I) {a)t - yp/3i 


C(t) = -2iSC{t) - {C{t) - mi) , 

v{t) = -iJL{u{t) - ni) . 

The solution to this set of equations reads 


{a)t = e-(*^+^)‘ 




pmi 


2i6 + fJ. J 2i6 + fi ’ 

iy{t) = {I'o — ni) + ni . 

Let us notice that the asymptotic solutions have the form: C(eo) = 2 ^^+^ ’ 
v{oo) = ni, and (a)oo = It means that the systems approaches 

asymptotically the Gaussian state with parameters independent of their initial 
values. It should be noted that for the resonance case {S = 0) the asymptotic a 
priori and a posteriori states coincide. 
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